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Asymptotics of Random Contractions 



ENKELEJD HASHORVAtl ANTHONY G. PAKES 0, and QIHE TANG£ 



^ ■ Abstract: In this paper we discuss the asymptotic behaviour of random contractions X = RS, 
^ . where R, with distribution function F, is a positive random variable independent of S G (0,1). 
! Random contractions appear naturally in insurance and finance. Our principal contribution is the 
derivation of the tail asymptotics of X assuming that F is in the max-domain of attraction of an 
extreme value distribution and the distribution function of S satisfies a regular variation property. 
p ( . We apply our result to derive the asymptotics of the probability of ruin for a particular discrete- 
^ ' time risk model. Further we quantify in our asymptotic setting the effect of the random scaling on 
£h the Conditional Tail Expectations, risk aggregation, and derive the joint asymptotic distribution of 
. linear combinations of random contractions. 
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1 Introduction 

Let R, S be two independent random variables with R > 0, S G (0,1) almost surely, and define 
X = RS. The random variable X is a random contraction of R via S. Random contractions or 
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random scalings are common in insurance and finance applications. Typically R models a random 
payment whereas S is a random discount factor. Several authors have studied random contractions 
in quite different contexts. Some recent contributions dealing with distributional and asymptotic 
properties of random contractions are Kotz and Nadarajah (2000), Galambos and Simonelli (2004), 
Gomes et al. (2004), Maulik and Resnick (2004), Tang and Tsitsiashvili (2003, 2004), Jessen and 
Mikosch (2006), D'Auria and Resnick (2006, 2008), Tang (2006, 2008), Denisov and Zwart (2007), 
Pakes and Navarro (2007), Resnick (2007), Beutner and Kamps (2008), Charpentier and Segers 
(2007, 2009), Hashorva (2008, 2009, 2010), Hashorva and Pakes (2010), Liu and Tang (2010). 

Our main goal in this paper is to investigate the tail asymptotics of random contractions when the 
tail asymptotics of R is known. An important motivation for this investigation is the fact that in 
insurance and finance applications assumptions often are made on the tail behaviour of a random 
payment modeled by R. If S represents the random discount factor applicable to the interval from 
the present to the payment time, then RS is the present value of the later payment R. In cases where 
the distribution function G of 5" is unknown, it is of some interest to know how the tail behaviour 
of the random contraction X = RS is determined by the corresponding asymptotic behaviours 
of the factors. One possible application is to approximating the Value at Risk in the presence of 
discounting, given information about the Value at Risk before discounting. 

Without going into mathematical details, we mention briefly the main contributions of this paper: 

a) Under the assumption that the distribution function F of R is in the max-domain of attraction 
of some univariate extreme value distribution we obtain the asymptotic behaviour of P{X > u} as 
u tends to the upper endpoint of F, provided that G = 1 — G satisfies a regular variation property 
(Theorem 13.11 below) . 

b) We determine corresponding results for the density function of X assuming a regular variation 
property for the density function of S, and additional regularity properties in some cases. 

c) We present four applications: cl) First we derive the asymptotics of the ruin probability for a 
particular discrete-time ruin problem, c2) then we discuss briefly the asymptotics of Conditional 
Tail Expectation in the random contraction framework, c3) and we obtain asymptotic expansions 
for the aggregation of two contractions, which lead to novel asymptotic characterisation of bivariate 
elliptical distributions, c4) finally we show the asymptotic independence of certain bivariate random 
contractions. 

As mentioned above, we assume that a generic scaling factor S with distribution function G takes 
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values in (0, 1). In addition, we assume that G(l— y) is regularly varying at zero, the above mentioned 
regular variation property. However, the reader will easily appreciate that the scaling factors can be 
multiplied by a positive constant and hence can function as inflation or deflation factors. Since our 
results can easily be adjusted for this contingency, we will say no more about it beyond the closure 
property in Lemma [2. II 

It is interesting that under the setup of this paper the asymptotic tail behaviours of R and X are very 
similar. In particular, membership of a max-domain of attraction is insensitive to the distribution 
of bounded discount factors. 

Our main results are presented in Section 3 followed by the applications in Sections 4. The proofs 
of all the results are relegated to Section 5. 



2 Maximal Domains of Attraction 

In this short section we present some details on max-domains of attraction. The distribution function 
F belongs to the max-domain of attraction of a univariate extreme value distribution function N, 
written F G MDA(iV), if 



lim sup 



F n (a n x + b n ) - N(x) = (2.T 



holds for some constants a n > 0, b n G M, n > 1. See e.g., Reiss (1989), Embrechts et al. (1997), 
Falk et al. (2004), De Haan and Ferreira (2006), or Resnick (2008) for more details on univariate 
max-domains of attraction. Only three choices for N are possible, namely the Frechet distribution, 
the Gumbel distribution, or the Weibull distribution. We denote the corresponding distribution 
functions by $ 7 , A, and \I/ 7 , respectively, where 7 > indexes members of the Frechet and Weibull 
families. 

The functional form of the Frechet distribution function is $ 7 (x) = exp(— x -7 ), x > 0. If F G 
MDA($ 7 ), then ([271]) with N = $ 7 is equivalent to 



l im = x -7 ; y x > 0. (2.2) 

u^oo F(u) 

This means that the survival function F = 1 — F is regularly varying at infinity with index —7 and 
further it has an infinite upper endpoint (denoted in the sequel by ?>). 
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The functional form of the standard Gumbel distribution function is A(x) = exp(— exp(— x)), x £ M, 
and (12. ip with N = A is equivalent to 

Um Fiu + xjw^ = yxeR ^ 

wtr F F(u) 

where w is a positive scaling function satisfying 

lim uw{u) = oo, and lim w{u){r F — u) = oo if r F < oo. (2.4) 

w\r F w\r F 

Recall that the scaling function w can be defined asymptotically via the mean excess function (see 
e.g., Embrechts et al. (1997) or Resnick (2008)) by 

mW ~ E{R-l\R>uY UUF - (2 ' 5) 
Throughout this paper he relation a(u) ~ b(u) means that the quotient of both sides tends to 1 
according to the indicated limit procedure. 

The functional form of the Weibull distribution function is ^ 7 (x) = exp(— |x| 7 ), x < 0. If F £ 
MDA(\l / 7 ), then r F is finite and (12.11) is equivalent to 

lim - X/U \ = x\ Vx > 0. (2.6) 

m^oo F(r F — 1/u) 

In some applications it is necessary to admit scaling factors S £ (0,c), where c is a positive con- 
stant. We always assume that c = 1. If F £ MDA(\l/ 7 ) we assume too that r F = 1. The next lemma 
explains why these conventions are not restrictive. 

Lemma 2.1. Let W be a random variable whose distribution function F satisfies (12.11) . If c,p £ 
(0, oo) , then cW p has a distribution function in the same max-domain of attraction as W . 

We will need the following facts about subexponential distribution functions. See Embrechts et al. 
(1997, Appendix A3) for distribution functions supported in [0, oo), and Borovkov and Borovkov 
(2008, p. 13) for the general case. Let F 2 * denote the convolution square of F, i.e., the distribution 
function of R + R*, where R* is an independent copy of R. Assume that r F = oo in what follows. 

In the case that F(0— ) = 0, we say that F is subexponential, written F £ S + , if 



lim = 2. (2.7) 

F(« 



u — ^oo 



In the case that F is two-sided, i.e., F(0—) > 0, define F + (u) = if u < and F + (u) = F{u) 
if u > 0. We say that F is subexponential if F + £ «S+, and we write F £ S; see Borovkov and 
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Borovkov (2008, p. 14). These authors (p. 19) show that it is possible that a two-sided distribution 
function F (jL S satisfies ( 12.71) . but imposing a condition described below will ensure that F G S is 
equivalent to ( 12. 1\ . 



Say that F belongs to the class of long-tailed distribution functions, written F G C, if 

Km ?t±i> = 1 

u^oo F(u) 

for all real y. The convergence here is locally uniform with respect to y. If F(0— ) = 0, (12. 7p implies 



that F G £. In the two-sided case, if F G £, then F 2 *{u) ~ F+*(u) as u — > oo (Borovkov and 
Borovkov (2008, Theorem 1.2.4(vi)), and hence F 6 5 is equivalent to (12. 7p . These concepts relate 
to attraction to the Gumbel distribution as follows. 

Still assuming that rp = oo, assume too that F G MDA(A, w) with lim u ^ 00 u;(n) = 0. For real y 
we can choose u so large that \y\ < \x\/w(u) and hence conclude from (I2.3P that F G C. Further 
conditions can be given to ensure that F G £ fl S. This holds if there is a positive constant A such 
that (see Mitra and Resnick (2008, Corollary 2.9)) 

Um [F(A/ mM F =a (28) 

Note that if F(0—) = 0, in view of Corollary 2.5 of Goldie and Resnick (1988), when w is eventually 
non-increasing such that 

lim 44 > 1 (2-9) 
holds for some constant £ > 1, then F G <S + . 

We derive by our next result a self-contained proof of the Mitra- Resnick criterion ( 12. 8p . 

Lemma 2.2. Let F G MDA(A, u;) with rp = oo and lim^oo w;(u) = (hence F £ £). Then F £ S 
if and only if 

pu—X/w(u) 

lim =— / F(u - y)dF{y) = (2.10) 

u^co F(u) J x/w(u) 

holds for some A > 0. 

By Lemma 12.21 the Mitra- Resnick criterion follows immediately, because the integral in (1 2 . 1 [) is 
bounded above by 

F(X/w(u)) \F(\/w(u))-F(u-\/w(u)j\ < [~F(\/w{u))} 2 . 
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Furthermore, Lemma 12.21 implies a generalization of Goldie's sufficient condition asserting that if 
F G C has a dominatedly varying right-hand tail, i.e., limsup u _ > . 00 F(u/2)/F(M) < oo, then F G S 
(see Embrechts et al. (1997, pp. 49, 52)). More precisely, suppose that F G MDA(A,w) with 
r F = oo and ]im u - ¥ . 00 w(u) = 0. Then by Lemma |2~2| F G S if 



= o (2.11; 



u— >oo 



for some A G (0, oo). 



3 Principal Results 

Let R be a positive random variable with distribution function F, and let Si, . . . ,S n be mutually 
independent, and independent of R. Denote by Gi the distribution function of Si, and assume it is 
supported on (0,1). In insurance and financial contexts Si represents the random discount factor 
over the interval [i — 1, i). Then 

n 

X n = R Y\ S i 

t=i 

is the present value of the payment R received at time n. Let H n denote the distribution function 
of the random product X n . 

If F G MDA($ 7 ) then, with no further conditions, it follows from Breiman's lemma (Breiman (1965)) 
that 

iim ^w4=Wri s -} (3-D 

and, in particular, that H n G MDA($ 7 ). See Jessen and Mikosch (2006), Denisov and Zwart (2007), 
and Resnick (2007) for details on Breiman's lemma and for some of its generalisations. So we need to 
consider only the cases where F is in the max-domain of attraction of the Gumbel or of the Weibull 
distribution. In these cases it turns out that the tail asymptotic behaviour of each Si is crucial. Our 
working assumption for the scaling random variables is that 

lim =x « % Vx>0 (3.2) 

«-*ooGi(l-l/«) 

for some aii G [0, oo). So if ctj > 0, then Gi G MDA( 1 l/ a J. If Gi possesses a positive density function 
gi, and if ctj G (0, oo), then we will assume the von Mises condition 

lim 9i ^~ X J, U l = x a >~\ Vx > 0. (3.3) 
u^oo gi (l - l/u) 
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The following is our first result, in which we denote by T(-) the Euler gamma function. 

Theorem 3.1. For i < n let Si (with distribution function Gi) be mutually independent scaling 
random variables, and independent of R (with distribution function F). Assume that F(0—) = 
and rp G (0, oo], and also that (13 .2D holds for every i <n with aij G [0, oo). 
a) If F G MDA(A, w) with r F G (0,oo], then 

n / 1 \ 

P{X n >u}~ \ \ Wipn + VjGi (1 —)F{u), utr F . (3.4) 

AJ-L V uw(u)J\ 



al) If, in addition, rp = oo, lim^oo w{u) = 0, and (12.101) holds, then H n G S + . 
b) If F G MDA(\& 7 ) with 7 G [0, 00) (hence rp = 1 fry onr convention), then 

P{X n >u}~ - r ^ w +1) f[ [r(«i + 1)G,( W )] ~F(u), u 1 1. 



(3.5) 



In the light of Theorem 13. 1[ if -F is in the Gumbel or the Weibull max-domain of attraction and 
each survival function Gi(l — u) is regularly varying at 0, then X n has a distribution function in 
the Gumbel or the Weibull max-domain of attraction, respectively. See also the proof of Theorem 
13.11 below. Part al) of Theorem 13.11 shows that the subexponentiality is preserved under random 
scaling. Recent results on the subexponentiality of products are obtained in Tang (2006, 2008), Liu 
and Tang (2010). 

The random contraction X n possesses a density function h n if one of the scaling random variables 
Si has a density function or if R has a density function /. In the following theorem we derive 
asymptotic approximations of the density function h n . Part a) is a density version of Breiman's 
lemma. The assumption under al) below places conditions on one of the density functions g^ but 
not on F (beyond A2.2D ) . In particular, can behave like a beta density near the origin, but it 
is bounded near unity. It seems that relaxing this condition requires restricting the slowly varying 
factor of F. One possibility is to assume that this factor is normalized slowly varying, i.e., that the 
density function / exists and is regularly varying. We do this in a2), and then we need no conditions 
on Gi. In parts b) and c) we assume for some i that g^ satisfies (13. 3 j) with an additional technical 
condition for b). 

Theorem 3.2. Let Si,i < n,n > 1, and R (with distribution function F) be as in Theorem \3.1\ 
a) If F satisfies (12. 2 p with 7 G [0,oo), then 

hm — — r- = 7 (3.6) 

«->oo P{X n > u} 
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holds under either of the following conditions: 

al) For some % <n, the scaling factor Si has a density function such that ygi{y) is bounded in (0, 1), 
or 

a2) F has a density function f which is regularly varying at infinity with index — (7 + 1) for some 
7 > 0, and Jq y~ e dGi(y) < 00 for some e > and all i < n provided that 7 = 0. 
In addition, lim^oo h n (u)/f(u) exists and equals the limit in (13. ip . 

b) Assume that f l3.2j) holds for all i < n with a, £ (0, 00), and F £ MDA(A, w) with rp £ (0, 00]. 
Suppose too for some i < n that Si has a density function g\ satisfying ( 13. 3 1) and that there exist 
c £ (0, 1) and Pi > such that 

sup y Pi gi(y) < 00. (3.7) 

0<j/<c 

Then 

h„(u) 

lim / NprV 7 = L 3.8 

c) Suppose that the distribution function F satisfies (I2.6P u>z£/i 7 £ [0, 00) (hence, rp = 1 by our 
convention) and that (13. 2ft /ioWs /or alii < n with oti > 0. // m addition ( 13.31) /ioWs /or some z < n 

nm — , r = 7+/ 3.9 

140 P{X„ > 1 - u} ' 4-f 

Under the assumptions of part b) of Theorem 13.21 we obtain further 

h n (u + x/w(u)) 

lim — — = exp(— x), Vx EJR. (3.10) 

v\r F h n (u) 

If F £ MDA(A, w) and it has a density function /, then conditions exist under which 

lim /(n + a:Mn)) =ex y 

see Resnick (2008). If ( 13. lip holds, then we can derive (13. lOf) under milder conditions than in The- 
orem [372](b) . In essence, we assume below that all the Oj exist and satisfy a condition similar to 
intermediate regular variation. 



Theorem 3.3. Let Si,i < n,n > 1, be mutually independent scaling random variables, with positive 
density functions gi, and independent of the random variable R > which has a distribution function 
F with an upper endpoint rp £ (0, 00]. Suppose further that F possesses a positive density function 
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/ such that (13.111) holds. If for i < n, s G (u,r F ), and any measurable function a : M 2 —> [0, oo) 
such that lim u p. F a(u, s) = 1 we have 



tc\r F gi[u/s) 



then f)3.10p is satisfied. 



Note in passing that if rp is finite and <fc(l — u) is regularly varying at with index on — 1 G (0, oo), 
then (I3.12p is satisfied. 

We present next three illustrating examples. 

Example 1. Let R be a random variable with distribution function F G MDA(A,u;) and upper 
endpoint rp G (0, oo], and let S be a random variable with beta distribution with positive parameters 
a, (3. Since 

P{g>l-«}^ F( °^ u", uiO, (3.13) 

it follows from Theorem 13.11 that the distribution function H of RS satisfies 

H{u) ~ r y HH)-^), u t r F . 

If g is the positive density function of S, then condition (13 .TP holds, whence (I3.8P implies that the 
density function h of if satisfies 

h(u) ~ w(u)H(u), u^rp. 

Note in passing that condition (I3.12p can be easily checked. 

Example 2. Under the setup of the previous example suppose that 

P{R > u] ~ Ku q exp(-ru~<), K > 0, r > 0, 7 > 0, q eM (3.14) 

holds for u — > 00. It follows easily that F G MDA(A,u>) with w(u) = r , ~fu' y ~ 1 ,u > 0. Hence relation 
(I3.13P implies 

H(u) ~ K(r 7 )~ Q ^"° 7 exp(-rn 7 ), u ^ 00. 
Further, by (I3.8p . the density function h of RS satisfies 

h{u) ~ r^u 1 ~ 1 P{RS > u}, u — > 00. 
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Note in passing that condition (\2.8\i does not hold for any A £ (0, oo). However, if 7 £ (0, 1), then 
R and RS have subexponential distributions. 

Example 3. Let R be a positive random variable with distribution function F satisfying 

F(u) ~ ciexp(-c 2 /(l -u)), ufl, 

with Ci, C2 two positive constants. Since for w(u) = c 2 /(l — u) 2 , u £ (0, 1), and any s £iR, 

F(u + s/w(u)) r ,. , , ... „ 

1 ' 1 " ~ exp(-c 2 [l/(l -u + s/w(u)) - 1/(1 - u) exp(-s), n | 1, 
F(u) 

we have F £ MDA(A,w). Let S* £ (0, 1) be a random variable independent of R such that (13.31) 
holds. Applying Theorem 13.11 we obtain 

H(u) ~ T(a + l)F(u)G ( 1 - (1 - ^ ] , ufl- 

4 Applications 

4.1 Ruin in the Presence of Risky Investments 

Consider the following discrete-time insurance risk model. Within period i, the insurer's net profit 
(total premium income less claim payment) is denoted by a real-valued random variable Z^. The 
insurer positions him/herself in a discrete-time financial market consisting of a risk-free bond with 
a constant periodic interest rate Si > and a risky stock with a periodic stochastic return rate A; 
taking values in (— l,oo). Suppose that, in the beginning of each period i, the insurer invests a 
fraction 7Tj £ [0, 1) of his current wealth in the stock and keeps the remaining wealth in the bond. 
Denote by Ui the insurer's wealth at time i, with a deterministic initial value Uq — u > 0. Then, Ui 
evolves according to 

Ui = [(1 - 7Ti)(l + Si) + 7T<(1 + A;)] + Z U 1= 1, 2, 



As usual, define the probability of ruin by time n as 



ib(u:ri) = P< min U < 



fo = « i n = l,2,.... 



Assume that Z\, Z 2 , ... are independent and identically distributed random variables, that Ai, A 2 , 
. . . are independent random variables, and that the two sequences {Zi, Z 2 , . . .} and {Ai, A 2 , . . .} 
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are mutually independent. Introduce 
1 



1 



1 + A; 



Z " Si ~ (1-^(1 + 50 + ^(1 + AO' Z_1 ' 2 ' 



(4.1) 



The random variable Tj is the random discount factor during period i of the risky asset and it takes 
values in (0, oo), the random variable Ri is the net loss during period i, and the random variable Si 
is the overall random discount factor during period % of the investment portfolio. Denote by F the 
common distribution of {Ri, i — 1, 2, . . .}. 

According to Tang and Vernic (2010), if F G S, then 



ijj{u;n) ~ ^pj_R fe Y[Si > uj, u^oo. (4.2) 

k=l i=l 

See related discussions in Tang and Tsitsiashvili (2003, 2004). The applicability of this formula 
requires explicit asymptotic expressions for the tail probabilities in f)4.2p and our main results clearly 
are crucial for this purpose. 

Notice from (14.11) that if P{Tj > it} is regularly varying at infinity with index —at for some 
<y.i G [0, oo), then P{Si > — l/u} is regularly varying at infinity with index — aii, where = 
(1 - 7r0 _1 (l + S^ 1 G (0, oo). Actually, 



»{s , i >S i -l/u}~p{T i >7r i a?u}, 



u — y oo. 



(4.3) 



Hence, if F G S D MDA(A, w) and for each 1 < i < n, the tail probability P{T» > u} is regularly 
varying at infinity with index — oti for some ccj G [0, oo), then applying Theorem 3.1(a) to relation 
(14.21) we obtain, with w fc = wlli=i V^*j 



1=1 



P\R k \ 

1 = 1 



^ («*) n 

i=l 
i=l L 



^ 

s, u k w(u k ) 



r (a, + 1 



— P{Tj > u fc w;(u fc )} 



M — ^ OO, 



where the last step is due to (14.3jl . We summarize all this as follows. 



Theorem 4.1. Consider the discrete-time risk model introduced above. If F G S fl MDA(A, iu) and 
/or eac/i 1 < i < n, the tail probability P{Tj > u} is regularly varying at infinity with index —a,i for 
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some cii G [0, oo), then 

n 



fc r 



fc=l t=l 



r (ag + 1; 



M — )■ OO. 



Note in passing that if F G MDA(A,w) with w such that lim^oo w(u) = 0, then in order to show 
that F G S we can utilise ( I2.10p . 



4.2 Asymptotics of Conditional Tail Expectation 

Let Si, i < n, R (with distribution function F), and X n be as in Theorem 13.11 If F G MDA(A, w), 
then to satisfies the self-neglecting property (see e.g., Reiss (1989) or Resnick (2008)) 

lim ™(" + *M")) = 1 (4.4) 

utr F w(u) 

uniformly with respect to z in every compact set of M. So if the conditions of Theorem 13.1( a) 
are satisfied, it follows from (13 .2p and (13.5P that H n G MDA(A, w). Consequently, we obtain the 
following asymptotic formula (recall ( 12. 5ft ) 

lim gffi " "]*■ - = 1. (4.5) 

In several insurance and finance applications the mean excess function is a crucial quantity (see 
Embrechts et al. (1997), p. 294). The result in ( 14. 5 j) shows that under the assumed conditions, the 
mean excess function is asymptotically invariant under random contractions. 

The Conditional Tail Expectation (CTE) for R with continuous distribution function is defined e.g., 
by 

CTE fl (u) := E{R\R > u] = E{R - u\R > u] + u, u> 0. 
Then in view of (14. 5p . 

lim CTE flM = 1 . m CTE x » =1+lim 1 =1 

u\r F U u^r F U u^rp UW(u) 

Consequently, if VaR p (X n ) denotes the Value at Risk (VaR) corresponding to the level p G (0, 1), 
i.e., 

VaR p (X n ) := inf {x : P{X n < x} > p} , 

then 

CTE Xn (VaR p (X n )) ~ VaR p (X n ), p f 1. 
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It is well-known that for continuous risks CTE is more conservative than VaR. The above asymptotics 
shows that in the Gumbel case CTE and VaR are asymptotically the same and that this relation is 
preserved under random scaling. 

4.3 Linear Combinations of Random Contractions 

In order to motivate the next applications we consider a bivariate scale mixture random vector 
(U\,U 2 ) with stochastic representation 

(U U U 2 ) = R{I x S,hVl^), (4.6) 

where R, with distribution function F, is almost surely positive, I\, I 2 assume values in {—1, 1}, and 
S, with distribution function G, is a scaling random variable taking values in (0, 1). Furthermore, 
suppose that Ii, I 2 , R, S are mutually independent. If S 2 follows a beta distribution with parameters 
1/2,1/2 and P{h = 1} = P{I 2 — 1} = 1/2, then (Ui, U 2 ) is a spherically distributed random 
vector, see Cambanis et al. (1981). (We shorten this by saying that (U±, U 2 ) is spherical). Hence by 
Lemma 6.1 of Berman (1983) 

Cl Ui + c 2 U 2 = ^c\ + c\XJ x , V Cl , c 2 ElR. (4.7) 

If S follows a beta distribution with parameters a, b, then (Ui, U 2 ) is a generalised Dirichlet random 
vector (see Hashorva et al. (2007)), and (14. 7p does not hold in general. 

Next, we derive the tail asymptotics of the aggregated risk 

U(p) = pUx + ^/l-p 2 U 2 , p G (0, 1) 

for some general scaling random variable S. If I\ = I 2 = 1, then U(p) is maximized with respect to 
S at S = p. Hence we make the following assumption about the local form of G at p: 

P{\S - p\ < t} = L p (t)t a » } a p G [0, 00), (4.8) 

for all t G (0, e), e > 0, where L p is positive and slowly varying at 0, and L p (0+) = if a p = 0. 
Clearly, if G possesses a density function g continuous at p, then ( 14. 8ft holds for any p G (0, 1) with 
a p = 1 and L p (t) = (2 + o(l))g(p) as 1 1 0. 
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Lemma 4.2. Let be two random variables taking values —1,1 with q = P{I\ = 1,I 2 — 

1} G (0, 1] an<i independent of a scaling random variable S with distribution function G. For given 
p G (0, 1) define a new random variable 



S(p) := phS + a/1 - p 2 I 2 VT^S*. 
If G satisfies (14. 8 j) u>rf/i L p (0+) = u>/ien a p = 0, i/ien 

> 1 - n} ~ g L p (v^)(2«(l - p 2 ))^/ 2 , u ; 0. (4.9) 

Note that if G is absolutely continuous with a positive density function g continuous at p, then we 
have 



P{S(p) > 1 - u} ~ ^(^8(1 - P 2 K u + (4.10) 

In view of Theorem 13. II and Lemma 14.2} we now study the tail behavior of U(p). First consider the 
Gumbel case i.e., F G MDA(A, «;) with rp G (0, oo]. Then with rj(u) = uw(u) for u > we have 

PMp) > u} ~ qT(a p /2 + l)^^^)" 1 / 2 ) f u t r*. (4.11) 



n(w) 

Since lim u ^ rF 7y(w) = oo it follows that lim u | rF P{U(p) > u}/F(u) = 0. If G possesses a density 
function g continuous at p, then 

P{f/(p) > ~ qT(l/2)g(p) ( ^T^ ) ^(u), u t r*. (4.12) 



?7(m 

Now if (Ui,U2) is spherical, then 04. Tp (obviously!) implies the tail equivalence 

P{U(p) > u} rsj P{U X > u}, u t r F . (4.13) 
The corresponding density function h p of U (p) satisfies 

h p {u) ~ ufr F . (4.14) 

The following converse result characterises spherical random vectors. 

Theorem 4.3. Suppose that Ii,I 2 are independent random variables assuming values —1,1 with 
probability 1/2, the distribution function F with F(0—) = satisfies ( 12. lft . and G possesses a 
continuous density function g. Then (I4.13P /io/ds /or a/Z p G (0, 1) i/ and only if (Ui, U 2 ) is spherical. 
Similarly, (14.141) holds for all p G (0, 1) if and only if (Ui, U 2 ) is spherical. 
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Next consider the Weibull case. Assume that F G MDA(\l/ 7 ) with 7 G (0, 00) (and rp = 1). 
Applying Theorem 13.11 and Lemma 14.21 we obtain 

P{U(p) >l-u}~ g r y 2 / ^ )r( ^/ ) L p (^)(2n(l - p^F(l -u), uiO. 

1 (pip/ 2 + 7 + I) 

If G possesses a continuous density function g, then this simplifies to 

p{u{ P ) > 1 - u} ~ /r^/a + y^ ^W 1 " W 2 ^ 1 - «)» u + °- 

The following result gives the Weibull analogue of Theorem 14.31 

Theorem 4.4. Suppose that I\,l2, g,G are as in Theorem \4-3\ If F G MDA(^ 7 ) with 7 G (0, 00) 
(and rp = 1), then (I4.13P holds for all p G (0, 1) if and only if (Ui, U2) is spherical. 

We remark that if F G MDA($ 7 ) with 7 G (0, 00), then the tail behaviour of U(p) follows from 
Breiman's lemma. Indeed, under this assumption U(p) has a distribution function in MDA(<3> 7 ). 



4.4 Max-Domain of Attraction of Bivariate Samples 

Suppose that Q is the distribution function of a bivariate random vector (X, Y). Extending ( 12. ID . we 
say that Q belongs to the max-domain of attraction of a bivariate max-stable distribution function 
iV if 



lim sup 



Q n (a n x + b n ,c n y + d n )-N(x,y) =0 (4.15) 



holds for some constants a n > 0, c n > 0, b n , d n G -K, n > 1. This implies that each univariate marginal 
distribution of Q is in the max-domain of attraction of the corresponding univariate marginal of N. 
Conversely, if further 

lim nP{X >b n ,Y > d n } = 0, 

n— >oo 

then (14.151) holds with N a product distribution function with univariate extreme value marginal 
distributions. 

In this last application we discuss the max-domain of attraction for the distribution function of 
(Ui, U(p)), defined in the previous subsection, assuming that F is in the max-domain of attraction of 
some univariate distribution. For any p G (0, 1) denote by Q p the distribution function of the bivari- 
ate random vector (Ui, U(p)), and let Qi lP , % — 1, 2, denote the corresponding marginal distribution 
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functions. We focus here on the cases where F is in either the Gumbel or the Weibull max-domain 
of attraction. Assume that the distribution function G of S satisfies (14.81) with a p G [0, oo), and L p 
positive, slowly varying at 0, and L p (0+) = if a p = 0. 

First consider the Gumbel case. Assume that F G MDA(A,u>). It follows from (14.111) that Qi iP G 
MDA(A, w) for % — 1,2. We show that, if p G [0, 1), then Q p is in the max-domain of attraction of a 
bivariate distribution function which is a product distribution. If rp is finite, then this follows from 
the evident fact that U\ and U (p) cannot simultaneously take their maximum values. If rp = oo, then 
set bi(n) := Q^ p (l — l/n),n > 1 with Q~j the generalized inverse of the i-th marginal distribution 
of Q p . It follows from Lemma 4.2 of Hashorva and Pakes (2010) that 

lim (&2( n ) — pbi(n)) = oo. (4.16) 

n— >oo 

Our claim follows if 

lim nPiUx > b^n), U(p) > b 2 (n)} = lim P{U(p) > 62(71) |Z7i > 6 x (n)} = 0. (4.17) 

n— >-oo n— >oo 

In view of Theorem 2 of Hashorva (2009) 

lim P\u(p) > p&i(ra)[l + x/y/b 1 {n)w{b 1 {n))]\U 1 > 61 (n)} G (0, 00), Vx eM 

implying (14. 17p . and hence our claim. Note that condition (I4.16P is satisfied for a distribution 
function F with tail asymptotics given by (13.141) . 

Next consider the Weibull case. Assume that F G MDA(\l/ 7 ) with 7 G (0, 00) and r> = 1. In view of 
fl4TTT]) it follows that Q hp G MDA(^ 7+ai/2 ) and Q 2 , p G MDA(^ 7+ap/2 ). Since r F = 1 it follows that 
Q p is in the max-domain of a bivariate distribution function which is a product distribution. By 
Lemma [2.11 this outcome can be formally generalized to the case where rp is an arbitrary positive 
constant. 



5 Proofs 



Proof of Lemma 12.11 In Lemma 5.2 of Hashorva and Pakes (2010) the proof of Lemma 12.11 is 
shown for F G MDA(A, w) and c = 1. The case c G (0, 00) and F is in the Frechet or Weibull 
max-domain of attraction can be easily shown, therefore omitted here. □ 
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Proof of Lemma [2^21 First of all, if relation (12. 10p holds for some A > 0, then it holds for all A > 0. 
This can be verified as follows. For all < Ai < A2 < 00, by F G MDA(A, w) we have, as u — > 00, 

"U— \i/w(u) 

F(u - y)dF(y) < F^/wiu)) {F(u - X 2 /w{u)) - F(u - \i/w{u)))] 

u—Xs/wfu) 

= o((e Aa -e Xl )F(u)) , 
and similarly, f£j™£) F ( u ~ v) dF (v) = { F ( u ))- 

Next, we assume F G C PI S. Since lim^oo w(u) = 0, it holds for every A > 0, every T > 0, and all 
large u > that 

pu—X/w(u) pT poo 

/ F(u - y)dF(y) < F^(u) - F{u- y)dF(y) - F{u- y)dF(y). 

J\/w(u) J~oo Ju+T 



By F G £ n S, it holds that F* 2 (u) ~ 2F(u), that 

f F(u-y)dF(y)~F(u)F(T), (5.1) 
J —00 



and that 

F(-T) < liminf =— / F(u - y)dF(y) < limsup =— / F(u - y)dF(y) < 1. (5.2) 

u^oo _P( M ) J u+T F(u) J u+ T 

Then, relation (I2.10p follows since T can be arbitrarily large. 

Finally, we assume that relation (12.101) holds for all A > and we prove F G S. This part of 
Lemma [2.21 is an easy consequence of Theorem 3.6 of Foss et al. (2009). Actually, for a distribution 
F G MDA(A,w;) with rp = 00 and lim^oo w{u) = 0, for any X u > with lim^oo X u = the 
function h(u) = X u /w(u) is F-insensitive in the sense of Foss et al. (2009), i.e., 

iim r(t L + *<"» = 1. 

Hence by Theorem 3.6 of Foss et al. (2009), relation (I2.10p implies F G S. Nevertheless, we give 
here another self-contained proof. Similarly as above, for every A > and T > we write 

(pT pX/w(u) pu—X/w(u) pu+T poo \ 

/ + / + / + / + / }F(u-y)dF(y). 

J -00 JT J X/w(u) Ju~X/w(u) Ju+T J 

Estimates of the first and last terms above have been given in (15 .11) and (15.21) . and an estimate of 
the third term is given by (12.100 . For the second and the fourth terms, we have, as u — > 00, 

1 f X/w(u) ^ , F(u-\/w(u)) 



F(u) 



[ U F(u~ y)dF{y) < F{U ^ {U)) (F(T) - F(\/w(u))) -+ e x F(T), 
Jt F(u) 
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and 



U+T w v 7 , * F(m — X/w(u)) — F(u + T) A 1 
F{u - y)dF(y) < ^ ^ ^ e A - 1. 



By the arbitrariness of A and T, we easily conclude that F* 2 (u) ~ 2F(u). □ 

Proof of Theorem [3J] a) The claim in (I3.4|) follows by Lemma A. 5 of Tang and Tsitsiashvili 
(2004) and Theorem 3 of Hashorva (2009). We give here another direct proof. Let X = RS, where 
the factors are independent and the distribution function G of S satisfies (13. 2p with the subscript i 
omitted, i.e., 

G(l -x) = x a L(l/x), (5.3) 
with a > and L slowly varying at infinity. Clearly, 

P{X >u}= I" G{u/y) dF{y). (5.4) 

J u 

Substitute y = u + z/w{u) and define the random variable W u by 

„ rTT , F(u + z/w(u)) „ . . 

P{W U >z}= K / K } \ < z < r(u), 

where r(u) = [rp — u)w(u) < oo if rp < oo, and r{u) = oo otherwise. Observe that (12. 3 j) can be 
expressed as W u -4 W,W\ r F , where W is a random variable following an exponential distribution 
with mean 1 and A means convergence in distribution. 

Next, (15 .4p can be expressed as 

P{X >u} = F(u)E{G(l/(l + W u /rj(u)))}, rj(u) = uw{u). 

The survival function in this expectation is asymptotically proportional to W£G(1 — l/r)(u)) almost 
surely as u t Tf- We consider two cases. 

i) If W u < 1, then W U /(W U + r)(u)) < 1/(1 + r)(u)), and hence dominated convergence implies that 

hm E{G(l-_W u / { W u + v{ u)) h W u <l} = w< 
«tr F G{1 - 1/(1 +7](u))) 

ii) If W u > 1, then Potter's bounds for slowly varying functions (Bingham et al. (1987, Theorem 
1.5.6 (i))) imply that for chosen constants A > 1 and 5 > 0, there exists a number r > such that, 
once r](u) > r (see (12.4|) ). 

C AW?.. 



L(l+ V (u)) 

Let M n denote the maximum of n independent copies of R. It follows from (12.11) that, for all 
positive k, the k-th order moments of \M n — b n \/a n converge to the fc-th order moment of the 
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Gumbel distribution (Pickands (1968)). This can be utilised to show that {W" +5 ;u < r F } is a 
uniformly integrable family. So it follows from (15. 3 p that 

E(Bil-W.nW. + «u)));W.>l} = w> 

^r F G(l-l/(l+77(«))) 

Since E{W a } = T(l + a), we conclude that, as u t rp, 

P{X > u} ~ T(l + a)G(l - l/r](u))F(u). 

The self- neglecting property of the scaling function (see (I4.4|) ) implies in addition that H e MDA(A, w) 
Hence H3.4[) can be proved by induction. 

al) The proof follows directly from Lemma [2.21 

b) Since r F = 1, we express (15 .4p as 

P{X>l-u}= f G{{l-u)/y)dF{y) = F{l-u)ElG[-^^)\ 

where 

P{W U < 2 } = ^j 1 "^) < ^ < 1. 
11 F(l-u) " " 

So as m —7- 0, we have W u A W 7 , where < z} = z 1 for < z < 1. Since 

1 - u x «(1 - W u ) 



1-«W U l-ulf„ ' 

it follows from (15 .3p that, almost surely, 

Clearly, for any u G (0, 1) 

b (^8 sB(i - m) ' 

hence dominated convergence yields that 



4oG(l-«) i \l-uWj) r(a + 7 + l) 



We prove (13. 5 p by induction as follows. Let 



and assume that 

"n-l 



#„_i(l - m) ~ C n _i 



n^(i-«) 



F(l - u). 
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Noting that the right-hand side is regularly varying at zero with index 7„_i, the case n—1 implies 
that 



# n ,(i- M )~a.„i£j(i-wvi); 

where P{VT (n _i) < z} = z^- 1 . But 



i=l 



F(l-u) 



{ J r («n + 7n-l + 1) 

and the assertion follows. □ 

Proof of Theorem 13.21 In all cases it suffices to prove the case n—1 because the general result 
follows from asymptotic estimates obtained from Theorem 13.11 applied to the product with n — 1 
contraction factors. Thus, omitting subscripts, the density function of X = RS is (see e.g., Lemma 
2.1 of Pakes and Navarro (2007)) 

h{u) = y~ 1 g(u/y)dF(y) = [ y' 1 f(u/y)dG(y), (5.5) 



o 



where the first equality applies when G has a density function, and the second equality applies when 
F has a density function. 



al) Let W u be a random variable whose distribution function is 

P{W u <z} = ^±, ze(o,l). 

Clearly W u A W as u — > oo, where P{W < z} = z 1 (so it is degenerate at provided that 7 = 0). 
Substituting y = u/z in the first integral of (I5.5P yields 

h{u) = ^M E {W u g(W u )}. 

u 

It follows from dominated convergence that the expectation converges to 

lim E{W u g(W u )} = E{Wg(W)} = 7 / g{z)z*dz = lE{S^}, 
and hence, from ( 13. II) . that 

uh(u] 

lim =^ = 7. (5.6) 

M-5.00 H(U) 

a2) The second integral form in ( 15.51) can be expressed as h(u) = -E{S' _1 /(mS'~ 1 )}. Expressing the 
regular variation assumption as f(u) = m~ 7 ~ 1 L(m), where L is slowly varying at infinity, we have 

^| = ,E;{^[L(n5- 1 )/L(u)]}. 
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Choose e G (0,7) if 7 > or choose as in the assumptions if 7 = 0. Since S < 1, it follows from 
Potter's bounds that for chosen A > 1 there exists u' > such that L(mS , ~ 1 )/L(m) < S~ € for all 
u > u' . Hence dominated convergence yields lim^oo h(u)/f(u) = E{S 1 }. The Karamata-Abelian 
theorem (Bingham et al. (1987, p. 26)) implies that ( 15. 6 j) still holds. 

b) As in the proof of Theorem 13.11 it suffices to consider the case n = 1, i.e., X = RS, where the 
density function g of S has the form 

g{l - l/u) = au- a+1 L(u), u > 1. (5.7) 

Letting rj(u) = uw(u), and recalling notation from the proof of Theorem 13. 1( a). it follows from ( 15. 5ft 
that the density function of X satisfies 

h(u) ~ - 1/t7(«))F(u)£?{(1 + WMu))- 1 !^}, u t r F , 

where 

= g<l - HV^yu) + w , a _ 1 / 1 + y(u) \ L(l + q(u)/Wu) 



(7(1-1/(7/(^ + 1)) " VW« + 77(U); L{l + Tj{u)) ' 

Note that (1 + W u /r](u)) < 1 and that the second and third factors of R u converge to unity almost 
surely. 

Let I be a (large) positive constant and B u = {W u < r)(u)/l}. Note that P{B U } — > 1 as u f rp. If 
a > 1, then the second factor of i? u is bounded above by (1 + l/rjlu)) ^ 1 —> 1, and if < a < 1, 
then on £? u the second factor equals 



W u + V {u)\ l - a ^ / (lti-'Mn) j < (l + rl) 



1 + r/(u) / \ 1 + 77(7/ 



1-0 



Next, choosing 5 G (0,max(a, 1)) and A > 1, I can be made so large that on B u and with u such 
that r?(w) > /, the third factor of R u is dominated by Potter's bound Amax(W^, W~ s ). Dominated 
convergence hence gives the conclusion 

lim E{(1 + Wjviu))-^; W u < rj(u)/l} = E{W a ~ 1 }. 

u-\r F 

Assume that rj(u) > I and consider outcomes on B u , i.e., that / < r)(u) < IW U . If a > 1, then 
the second factor of _R U is bounded above by unity, and if < a < 1 then this factor is dominated 
by W^~ a . To deal with the third factor, observe that condition ( 13.71) is equivalent to the existence 
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of p > such that L(l + z) — 0(z~ p ) as z — > 0. Consequently, I can be chosen so large that the 
numerator in the third factor of R u is O p [(W u /r](u)) p ], as « | tf- It follows from these estimates 
that 

(FJT/T/ p+max ( a " 1,0 )- R \\ 
fa('")) p ^fa('")) J ~ > °' «tr, 

since the denominator tends to infinity. 

Since EiW"- 1 } = T(a), it follows that, 

h{u) ~ Y(a)u~ 1 g{\ - l/r)(u))F{u). 

But g(l — y) ~ (a / y)G(l — y) , so (13. Sp (with the subscripts omitted) follows. The general result 
follows from Theorem 13 . 1 1 after replacing F with the distribution function of RYljj*, Sj- 

c) The proof for the case n — 1 is similar to those above. Observe that if u,w £ (0,1), then 
(1 — uw)/ (1 — w) > 1. Assuming (15. 7p . then choosing 5 £ (0, a), we have Potter's bound 

w - - ^») £ ^ max(lj (1 _ uwy4(1 _ WnY% u > u , 

where W u is as in the proof of Theorem 13.1( b). This leads to the asymptotic form 

h{\ - u) ~ g (l - U )F(1 - u)E{(l - W) a ~ 1 }, u | 0. 
The assertion for n = 1 then follows with no further assumptions about g. 

For the general case we may assume with no loss of generality that (15. 7p holds with i = n, and then 
apply the single factor case with H n _x replacing F. It follows from Theorem 13.1( b) that 

h n (l - u) ~ g n (l - u)ff n _i(l - u)E{(l - Wn-x) *- 1 }. 

But g n (l — u) ~ (a n /u)G n {\ — u) and 

„ r . jpjn ttt \a n -x x n «nr(Q w )r(7 ra _i + 1) 

r (7n) 

It follows that /i n (l — u) ~ 7„if n (l — tt)/w, whence the general result (13. 9p . □ 

Proof of Theorem 13.31 It suffices to consider the case n = 1. We can, with no loss in generality, 
choose the scaling function to be different iable and satisfy 

_. . d 1 

u;(s) := -> 0, utr>. 

as w s 



23 



Replacing u in (15.5P with u + x/w(u), the substitution y = s + x/w(s) yields 

,/ / / w C F fu + x/w(u)\ 1+W(s) i / \\ j 

h(u + x w(u)) = gl ' )[ \\ f(s + x w(s))ds. 

J u \s + x/w(s) J s + x/w(s) 

Write the argument of g as (u/s)a(u, s), where 

. . l + x/(uw(u)) . 
1 + xj {sw{s)) 

Also, the middle factor in the above integrand is asymptotically proportional to Hence (I3.10P 
follows from ( 13.1 ip and ( I3.12P ; see Lemma 5.1 of Takahashi and Sibuya (1998). □ 

Proof of Lemma 14.21 Some algebra shows that S(p) < 1 and it is bounded away from unity unless 
I\ — I2 — 1. If this event occurs, then S(p) is close to 1 if and only if S is close to p. So on the 
event {I\ — I 2 — 1}, algebra reveals that, for all small u > 0, S(p) > 1 — u if and only if 

(p - S) 2 + 2pSu <2u- u 2 , 

and this condition is equivalent to 

(p - S) 2 < (1 + o p (l))2(l - p 2 )u. 

But 

P{{p ~ Sf < t} ~ t a ^ 2 L p (Vi), u I 
and the assertion follows. □ 

Proof of Theorem S3] The assumptions imply that (14.121) holds for all p e (0,1). The tail 
equivalence ( I4.13P implies that the factor g(p)(l — p 2 ) 1 ^ 2 is constant with respect to p G (0, 1). The 
density function of S 2 is g(y 1 ^ 2 ) /2y 1 ^ 2 , i.e., S 2 has the beta distribution with parameters (1/2, 1/2). 
If ( I4.14p holds, then so does ( 14 . 1 3 j) . Hence either condition implies that (Ui, U2) is spherical. □ 

Proof of Theorem S3] The proof is similar to that of Theorem S23 a 
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